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We have analyzed the NA58 data on J/t/j suppression in Pb+Pb colhsions. J/ip production is 
. . . assumed to be a two step process, (i) formation of cc pair, which is accurately calculable in QCD 

' and (ii) formation of J/tp meson from the cc pair, which can be conveniently parameterized. In a 

^— I pA/AA collision, a cc pair gain relative square momentum as it passes through the nuclear medium 

and some of the cc pairs can gain enough momentum to cross the threshold to become an open charm 
meson, leading to suppression in pA/AA collisions. A new prescription is proposed for the gain in 



o 
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^ , momentum square, consistent with Krammer process. The model without any free parameter could 

■ explain the Et dependence of J/ip over Drell-Yan ratio. 

oo . 

In relativistic heavy ion collisions J/ip suppression has been recognized as an important tool to identify the possible 
phase transition to quark-gluon plasma. Because of the large mass of the charm quarks, cc pairs are produced on 
• a short time scale. Their tight binding also makes them immune to final state interactions. Their evolution probes 
' the state of matter in the early stage of the collisions. Matsui and Satz predicted that in presence of quark-gluon 
. plasma (QGP), binding of cc pairs into a, J/ip meson will be hindered, leading to the so called J/ip suppression in 
' heavy ion collisions Over the years several groups have measured the J/ip yield in heavy ion collisions (for a 
, review of the data and the interpretations see Refs. ||^,^). In brief, experimental data do show suppression, which 
' could be attributed to the conventional nuclear absorption, also present in pA collisions. 

The latest data obtained by the NA50 collaboration Q on J/^ production in Pb+Pb collisions at 158 A GeV is 
■~>«^ ■ the first indication of the anomalous mechanism of charmonium suppression, which goes beyond the conventional 
I suppression in a nuclear environment. The ratio of J/ip yield to that of Drell-Yan pairs decreases faster with Et 

pattern can be 
Essentially 
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I . in the most central collisions than in the less central ones. It has been suggested that the resulting 
^ understood in a deconfinement scenario in terms of successive melting of charmonium bound states 
^ , in a QGP like scenari, assuming all the Jhp melts above a threshold density, NA50 data could be explained as an 
^ ■ effect of transverse energy fluctuations . The data could be also be explained in the comover approach without 
invoking QGP like scenario Recently we have shown that the NA50 data could be well explained extending the 
model of Qiu et al [lO| ]. In this model suppression is due to gain in relative square momentum of cc pair, as it pass 
through the nuclear medium. Some of the cc pair might gain enough momentum to cross the open charm threshold. 
. In the present paper, we further analyze the problem and alter the prescription used by Qiu et al [|o| to calculate 
" " ' the gain in the relative square momentum of the cc pair to be consistent with theoretical understanding. We argue 
that square of relative momentum gained by cc pair for traversing a length L in nuclear medium goes as rather 
than L, as was proposed by Qiu et al [0. 

We briefly describe the model of Qiu et al ||l^. Production of J/ip meson is assumed to be a two step process, (i) 
production of cc pairs with relative momentum square g^, and (ii) formation of J/ip mesons from the cc pairs. Step 
(i) can be accurately calculated in QCD [ p^[pd| ]. The second step, formation of J/ip mesons from initially compact cc 
pairs is non-perturbative. They used a parametric form for the step (ii), formation of J/ip from cc pairs. The J/ip 
cross section in AB collisions, at centre of mass energy y/s was then written as, 

CTA+B^J/4,+xis) = J '^1^ (-^^T^j J dxF(t)a/AiXa:Q'^)(l>b/B{xb,Q'^)^r^^^ X Fcc^j/^{q'^), (1) 

a.b 

where J2a b runs over all parton flavors, and = + 4m^. The K factor takes into account the higher order 
corrections. The incoming parton momentum fractions are fixed by kinematics and are Xa — i\/ x'p + 'iQ'^/s + xf)/'2, 
and xt, = {\Jx\ + / s — xp)/2. Subprocess cross sections can be found in ||l|] 

In a nucleon- nucleus/nucleus- nucleus collision, the produced cc pairs interact with nuclear medium before they exit. 
Observed anomalous nuclear enhancement of the momentum imbalance in dijet production led Qiu, Vary and Zhang 
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JlOl to argue that the interaction of a cc pair with nuclear environment, increases the square of the relative momentum 
between the cc pair. They prescribed that if the pair traverses a length L in nuclear medium, then the relative square 
momentum in the transition probability Fcc^j/^{q^) should be changed to, 

^q^+ e^L (2) 

with being the square of relative momentum gained by the pair per unit length of nuclear medium. However, if we 
treat momentum gain as a random walk, then for a Krammer like process, one easily obtain, 

<{pf-p^f >(x<{xf-x^)^ > (3) 

In other words, for a Krammer like process, square of momentum gain goes as rather than L as in Eqj^. 
Accordingly, we suggest that for traversing a length L in nuclear medium, the relative momentum q^ in the transition 
probability Ff.c^j/^{q^) should be changed to 

q^^q^+e^L^ (4) 

The parameter was obtained by fitting NA50 data on J/ifj production in pA and AA collisions. For the 
transition probability, we use the following form [p^,p[. 



Fcc^j/^{q^) = Nj/^e{q^)e{Am'^ - Ami ~ <f){^ - , / , J "", (5) 

Am' ~ Ami 



and we have used CTEQ5M parton distribution function in the calculation. In Fig.l, we have compared NA50 
experimental data with the best fitted curve. Except for the pd reaction, all the data points are correctly 
reproduced. Cross section for pd reaction is underpredicted. For a Krammer process, EqJ| (Ap^ c>c Aa;^) is obtained 
with the assumption of long interval of time i.e. for large interval of length. For pd collisions, effective length is small 
for the equation to be valid. 

We now apply the model to obtain transverse energy dependence of the J /ip over Drell-Yan ratio. The Drell-Yan 
pairs do not suffer final state interactions and the cross section at an impact parameter b as a function of Et can be 
written as, 

(fa'^'' /dETd% = a^l J d^sTA{s)TB{s - h)P{b, Et), (6) 

where cr^Jj is the Drell-Yan cross section in NN collisions. All the nuclear information is contained in the nuclear 
thickness function, TA,B{^)i= J dzpA.sis, z). Presently we have used the following parametric form for pa{^) @, 

PA{r) = C^) 

with a = 0.53fm, tq = lAA^^^. The central density is obtained from J pA{r)d^r = A. In Eq.^, P{b,ET) is the 
probability to obtain Ex at an impact parameter b. Geometric model has been quite successful in explaining the 
transverse energy as well as multiplicity distributions in AA collisions Transverse energy distribution in 

Pb+Pb collisions also could be described in this model In this model, Et distribution is written in terms of 
Et distribution in NN collisions. One also assume that the Gamma distribution, with parameters a and (3 describe 
the Et distributions in NN collisions. Pb+Pb data on Et distribution could be fitted with a = 3.46 ± 0.19 and 
/? = 0.379 ±0.021 |l. 

While Drell-Yan pairs do not suffer interactions with nuclear matter, the J/tp mesons do. In the model suppression 
factor depend on the length traversed by the cc mesons in nuclear medium. Consequently, we write the J/tp cross 
section at an impact parameter b as, 

d^a-^^^/dETd^b = ail^^ J d^sTA{s)TB{s - h)S{L{h, s))P{b, Et), (8) 

where cr'l^^ is the J/ip cross section in NN collisions and S{L{h, s)) is the suppression factor due to passage through 
a length L in nuclear environment. At an impact parameter b and at point s, the transverse density can be calculated 

as, 

n(b, s) = Ta{s)[1 - e-'^""'^«(b-s)] _^ j,^^^ _ ^^j^ _ g-<T„„TA(s)]^ (g) 
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and the length L(b, s) that the J/ip meson will traverse can be obtained as, 

L(b,s) =n(b,s)/2po (10) 
Suppression factor S[L{h,s) is then easily calculated using Eqj^. 

Fluctuations of Et, at a fixed impact parameter plays an important role in the explanation of NA50 data. Following 
Blaizot et al we take into account Et fluctuations by the following replacement, 

L{h,s)^ L{h,s)ETlET{b) (11) 

In Fig. 2, we have compared the model prediction with the NA50 experimental data on J /ip over DY ratio. The 
solid line is the ratio obtained with out any Et fluctuations. The data are reproduced upto the knee of the Et 
distribution. Beyond the knee of the Et distribution, model prediction saturates, while the data shows rapid fall of 
the ratio. The dashes line is the model prediction obtained including the Et fluctuations. The model without any free 
parameter give excellent description of the data, throughout the Et range. The results clearly shows that the rapid 
fall of the J/-0 over Drell-Yan ratio beyond the knee of the Et distribution is due to Et fluctuation only. Present 
calculation also establishes that it is not essential to assumed a deconflned scenario to explain the NA50 data on J /'ip 
to Drell-Yan ratio. Nuclear absorption alone is capable of explaining it. 

To summarize, we have analyzed the NA50 data on transverse energy distribution of J /ip in Pb+Pb collisions at 
CERN SPS. In the model, cc pairs gain relative square momentum as it travel a square of length through the 
nuclear environment. Suppression occurs as some of the pairs can gain enough momentum to cross the threshold to 
become an open charm meson. The parameters of the model were fixed by fitting experimental J /ip cross section in 
pA and AA collisions. The model could very well describe the transverse energy dependence of J /-ip over Drell-Yan 
ratio upto the knee of the Et distribution. It explains the data, throughout the Et range if Et fiuctuations are 
included. 
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FIG. 1. Total J/i/) cross sections with the branching ratio to H'^H in proton-nucleus, proton-nucleus and nucleus-nucleus 
collisions, as a function of the effective nuclear length L(A,B). 




FIG. 2. Open and closed circles are the J/ip to Drell-Yan ratio in a Pb+Pb collision obtained by NA50 collaboration in 
1996 and 1998 respectively. The solid line is the model prediction obtained without any Et fluctuations and the dashed line is 
obtained including Et fluctuations. 
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